Multiple Scattering of Electromagnetic Waves by an Array of Parallel Gyrotropic Rods 
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We study multiple scattering of electromagnetic waves by an array of parallel gyrotropic circular rods and 
show that such an array can exhibit fairly unusual scattering properties and provide, under certain conditions, 
a giant enhancement of the scattered field. Among the scattering patterns of such an array at its resonant 
frequencies, the most amazing is the distribution of the total field in the form of a perfect self- similar structure 
of chessboard type. The scattering characteristics of the array are found to be essentially determined by the 
resonant properties of its gyrotropic elements and cannot be realized for arrays of nongyrotropic rods. It is 
expected that the results obtained can lead to a wide variety of practical applications. 
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Multiple scattering of electromagnetic waves by period- 
ically spaced elements demonstrates many intriguing fea- 
tures that are of considerable practical and scientific signifi- 
cance OJE). Knowledge of the individual properties and mu- 
tual location of elements in a periodic structure illuminated 
by an incident electromagnetic wave makes it possible to de- 
termine the spatial distribution and other characteristics of the 
scattered field. At a fixed frequency, the properties of periodic 
structures can be controlled by a change in the positions and 
material parameters of the scattering elements. In many cases, 
it is almost impossible to change the internal geometry of a pe- 
riodic structure such as a photonic crystal, for example, with- 
out its destruction or serious damage. However, the properties 
of the scattering elements in the array can easily be controlled 
if they consist of resonant gyrotropic materials, the parame- 
ters of which in some frequency ranges can be very sensitive 
to even slight variations in an external dc magnetic field. Such 
a possibility opens up new promising prospects for control- 
ling the scattering characteristics of an array of such elements 
without changing their sizes or positions. Despite significant 
progress in the analysis of multiple scattering by periodically 
spaced resonant nongyrotropic elements (6||3, the results ob- 
tained for them cannot be used for arrays of gyrotropic scatter- 
ers because of the fundamental difference between the scatter- 
ing characteristics of nongyrotropic and gyrotropic elements. 
Although gyrotropic photonic crystals have recently been dis- 
cussed in the literature OH), most works on the subject do 
no consider the role of individual resonant properties of the 
periodically spaced gyrotropic elements constituting the array 
in the formation of the scattering pattern. Meanwhile, one can 
expect that the joint contribution of the individual and collec- 
tive resonance scattering mechanisms to the diffracted field 
of arrays containing gyrotropic elements should lead to many 
interesting phenomena, the features of which are yet to be de- 
termined. In this work, the scattering of a normally incident 
plane electromagnetic wave by a two-dimensional array con- 
sisting of parallel resonant gyrotropic circular rods is consid- 
ered and the unusual scattering properties of such a structure 
are revealed and discussed. 

Consider an equidistant array of identical parallel rods of 



radius a (see Fig. [T]). The rods are embedded in a uniform 
background medium and aligned with an external magnetic 
field which is parallel to the z-axis of a Cartesian coordi- 
nate system (x,y,z). The axes of the rods lie in the xz 
plane and are specified by the relations x = jL and y = 0, 
where L > 2a and j = 0, ±1, ±2, .... The medium in- 
side each rod is described by the permittivity tensor e which 
is typical of a magnetoplasma and has the following nonzero 
elements: e pp = = e e, e p(f) = -e^ p = -ie g, 

and e zz = e^r]. Here, eo is the permittivity of free space, 
e = l- uj 2 p /(uj 2 - g = uIu h /[(uj 2 - and 
r] = 1 — u 2 /u 2 , where uo p and u H are the plasma frequency 
and the gyrofrequency of electrons, respectively, and u is the 
angular frequency. The medium outside the rods is isotropic 
and has the dielectric permittivity e out = e^e. The incident 
wave is assumed to be a TE monochromatic plane wave whose 
magnetic field is polarized in the z-direction. We will not con- 
sider the incidence of a TM wave because its scattering is 
unaffected by the gyrotropic properties of the rods and sim- 
ilar to that in the case of isotropic rods. The wave vector 
k of the incident wave has the components k x = —kcosO, 
k y = — ksinO, and k z = 0, where k = k^e 1 ! 2 is the wave 
number in the surrounding medium (ko is the wave number in 
free space) and is the angle of incidence, which is reckoned 
from the positive direction of the x-axis. 

Omitting a time factor of exp(icjt), the total magnetic field 
normalized to the incident-wave amplitude can be sought in 
the form 

OO CO 

H z = e~ ikr + J2 E Dj, m H^(k Pj )e- im ^, (1) 

j = — 00 171— — OO 

where m is the azimuthal index (m = 0, ±1, ±2, ...), Dj^ m is 
the multipole coefficient characterizing the scattering by the 

(2) 

jth rod to the mth azimuthal harmonic of the field, Hm is 
the Hankel function of the second kind of order m, pj is the 
distance from the axis of the jth rod to the observation point 
in the incidence plane, and <j)j = arcsin(?//p J ). In Eq. |T|), the 
first term represents the field of the incident wave, while the 
other terms account for the scattered field. 



2 



y 



2a 




©- 



V 



Figure 1 . (color online) Geometry of problem. 



The magnetic field H Z] j i7n of the mth azimuthal harmonic 



inside the jth rod is written as 



H Z ;j,m — Bj^ rn J rn {hQqpj^j 1 



(2) 
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where S m is the single-rod scattering coefficient: 

J^(Q)Jm(Q)-e^ 2 J m (Q)Srr 



H^(Q)J m (Q) - &I 2 H%\Q)£ S 
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Here, 6 m = (e 2 - g 2 )' 1 [eqJ^(Q) + mg(k Q a)- 1 J m (Q)} , 
Q = koqa, Q = hoe 1 / 2 a, and the prime denotes the derivative 
with respect to the argument. 

The translational symmetry of the problem makes it possi- 
ble to use the discrete Fourier transform (with respect to j) 
and its inverse: 
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Dj, m = D m (h)e ihL *dh. 
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(5) 



The Fourier transform of the incident- wave field comprises 
the Dirac function 5(hL — kL cos 6). Therefore, we will seek 
D m {h) in the form D m (h) = D m 5(hL — kL cos 6). As a 
result, we arrive at the following system of equations for D m : 
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where 

CO 

G m =^H$(kLl)[e ikLlc ° sG + (-l)™ e -ikLlcos6 
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If the rods are electrically small such that h <C 1, we can 
restrict ourselves to the dipole approximation and retain only 
the terms with m = ±1 and n = ±1 in Eq. ([6]). Then from 
Eq. ^ we have 



±i 



[(S'l - G 2 )e^ e - Goer^} 
(SZI-G 2 )(S^-G 2 )-G 2 



(8) 



where the multipole coefficient B^ m characterizes the field of ^ 
the corresponding harmonic, q = [{e 2 — g 2 )/s] 1 ^ 2 , and J m is 
the Bessel function of the first kind of order m. 

Satisfying the boundary conditions for the tangential field 
components H z .j, m and ; j, m on the surface of the jth rod 
and using the standard technique based on the scattering ma- 
trix method (HE), we can exclude the coefficients Bj j7n and 
obtain a system of equations for the scattering coefficients 
Dj^ m in the form 
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Figure 2. (color online) Frequency dependences of the real and imag- 
inary parts (solid and dashed curves, respectively) of the array factors 
(a) Go and (b) G2, the single-rod scattering coefficients (c) S-i and 
(d) Si, and the coefficients (e) D-\ and (f) D\. The resonances of 
Co, 2 are labeled in order of increasing frequency. The resonances of 
S-i and Si at the frequencies uji and cj u and the corresponding reso- 
nances of D-i and D\ are denoted as I and II, respectively. The val- 
ues of parameters are chosen to be Up/un — 6.47, u p a/c = 0.18, 
L/a = 132.8, e = 1, and = tt/2. 

The array factors Go and G 2 in Eq. ([8| have resonances in 
the cases kL cos = 2ttv for |0| < tt/2 and kL — 2iru for 
\6\ = 7r/2, where v = 1,2,... . We will denote the corre- 
sponding resonant frequencies by uo v . Figure [51 shows the fre- 
quency dependences of the quantities in Eq. for = tt/2. 
It is seen in Figs.[2ja) and[2jb) that Go and G 2 are discontin- 
ues functions of kL. When approaching the resonant value of 
kL from the side of larger kL, the real and imaginary parts 
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of G m tend to infinity and certain finite values, respectively. 
But if the resonant value of kL is approached from the side of 
smaller kL, then the real and imaginary parts of G m tend to 
some finite values and infinity, respectively. Such behavior no 
longer takes place if the array is nonequidistant, the number of 
rods in the array is finite, or the surrounding medium is lossy. 

In the case where k <C 1 and ko\q\a <C 1, the single- 
rod scattering coefficients S-i and Si have the resonant 
frequencies uj 1 and lj u , respectively, which are located on 
different sides of the surface plasmon resonant frequency 
uj r of an isotropic plasma column [9]. The frequencies 
uj 1 and oj u are approximately determined by the equation 
( £ +9Xt) [l + (^o&) 2 /2] = —e, where x T = ±1 for rn = 
=F±. It can be shown from this equation that under the addi- 
tional condition e = 1, one obtains uj r = uj p /V2, uj p /2 < 
uo Y < (jj p /y/2, and lj u > uj p /V2. With decreasing external 
magnetic field (uj h — > 0), the resonant frequencies uj 1 and uj u 
tend to uj p /y/2. With increasing external magnetic field (for 
uj h ^> uj p ), uj 1 — » ujp/2 and uj u — » uj h . 

Interaction of the individual and collective mechanisms of 
scattering can lead to both an increase and decrease in the scat- 
tering from the array compared with the scattering by a single 
rod. To demonstrate this fact, Figs.[2je) and[2jf) show the fre- 
quency dependences of D±i in the case where the frequency 
cos of the v = 3 collective resonance of the array coincides 
with the resonant frequency uj u of a single gyrotropic rod. In 
this case, the coefficient D_i has the most pronounced res- 
onance peak at a frequency that is very close to, but slightly 
lower than the frequency lj u along with much less pronounced 
peaks at oj u and uj v . At the same time, the coefficient £>i has 
minor peaks near uj 1 and uj v . However, the behavior of Di in 
the vicinity of uj u differs significantly from that of D-i near 
uo l . An important feature of the coefficients D±i is that they 
decrease and become comparable, so that the scattering occurs 
as in the case of isotropic rods, if the frequency uj tends to any 
of the quantities uj v from the side of the higher frequencies. In 
this limit, the coefficients D±i turn out to be very small and 
the array becomes almost transparent for the incident radia- 
tion. 

Generally, the field scattered by an array of gyrotropic rods 
differs significantly from that in the case where an array con- 
sist of isotropic cylindrical scatterers. Before proceeding in 
comparison of these two cases, we turn to the fields scattered 
by a single rod when the incidence angle 6 = it/ 2. Fig- 
ures |3ja) and [3jb) show the snapshots of the scattered-field 
component H s z c at the frequencies uj 1 and uj u in the case of 
a single gyrotropic rod. Each of the presented fields has a 
pronounced helical structure and differs significantly from the 
field scattered by an isotropic rod at the resonant frequency 
uj = uj y [see Fig.[3jc)]. It is seen that an isotropic rod scatters 
the incident wave rather weakly in the direction perpendicular 
to the wave vector of the incident wave, whereas the scatter- 
ing pattern of a gyrotropic rod is much more uniform. Far 
from the rod, this pattern resembles the spatial structure of 
the m = azimuthal harmonic of the field scattered by an 



isotropic rod [see Fig.[3jd)]. 




Figure 3. (color online) Snapshots of the scattered field H s z c in the 
cases where a plane waves is incident at an angle — tt/2 on a 
single gyrotropic rod (j = 0) at the frequencies (a) uj = ujj ~ 4cj h 
and (b) uj = uju ~ 5cj h and (c) on an isotropic rod at the frequency 
uj = u T . (d) The snapshot of the m = azimuthal harmonic of the 
field scattered by an isotropic rod at uj — uj t . The parameters cj p /cj h , 
ujpd/c, and e are the same as in Fig. 2. 

The above-mentioned features directly affect the multiple 
scattering by an array of cylindrical rods. For example, the ab- 
solute value of the total field H z shown in Fig.[4ja) for the case 
of scattering of a plane wave at the frequency uj = uj u ^ uj^ 
by the array of gyrotropic rods has a periodic structure in 
the form of a chessboard. The absolute value of the total 
field in the case of scattering by the array of isotropic rods 
at cj = uj r — CJ5 does not possess such a spatial structure [see 
Fig. |4jb)]. To better clarify the difference between the both 
cases, it is instructive to compare the corresponding scattered 
fields, which are shown in Fig. |4jc) and |4jd). The absolute 
value of the field scattered by the array of gyrotropic rods is 
periodic along the x-axis and slowly varies in the direction of 
the y-axis. This is related to an almost axisymmetric far-field 
scattering pattern of each of the gyrotropic rods forming such 
an array. If the field scattered by an array of isotropic rods 
were dominated by the m = azimuthal harmonic, then this 
field would by very similar to that shown in Fig.Qc). How- 
ever, such a spatial structure cannot by observed for the array 
consisting of isotropic rods, because their individual scattered 
fields are of dipole nature and form another pattern presented 
inFig.gd). 

Under certain conditions, one of the coefficients D-i and 
Di can increase extremely greatly, so that the total field is pre- 
dominantly determined by the scattered field. This situation 
is depicted in Fig. [5] and takes place when both the real and 
imaginary parts of the denominator in Eq. ^ simultaneously 
tend to zero. Such a giant enhancement of the scattered field 
is observed in the case of a small mismatch between any of 
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Figure 4. (color online) Absolute value of the total magnetic field 
H z for the arrays of (a) gyrotropic and (b) isotropic rods at uj = uj u 
and uj = (jj r , respectively, if kL/2n = 4.999, L/a = 219.565, and 
— 7r/2. The absolute value of the scattered magnetic field H s z c 
for the arrays of (c) gyrotropic and (d) isotropic rods at the respec- 
tive frequencies indicated above. The same values of the parameters 
uj p /uJh, uj p cl/c, and e as in Fig. [5] 



the resonant frequencies uj v of the array and one of the single- 
rod resonant frequencies (uj 1 or uj u ), provided that uj is close 
to the corresponding frequency uj v . For the plot of Fig. [5] 

uj — ujj — CJ4 and \D-i\ ^> 
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Figure 6. (color online) Snapshots of the scattered field H s z c in the 
cases of oblique incidence for (a) = 3tv/7 and (b) = 7r/7, when 
kL/2-K = 1.0002. The values of other parameters are the same as in 
Fig.|2] 



monic in the case of wave scattering by isotropic periodic 
structures ifTOlfTTIl . 

In conclusion, we note that the results obtained can be use- 
ful in creating media with 2D distributed feedback fT2l . de- 
veloping promising photolithography methods 1 13 ], operating 
multi-tube helicon plasma sources fT4lL and understanding the 
features of wave scattering by periodic magnetic-field-aligned 
plasma density irregularities in the ionosphere |[T5l[T6l . 
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Figure 5. (color online) Snapshots of the total magnetic field H z 
for the array of gyrotropic rods at uj = 0.998cji (uji = 4.0645cj h ), 
kL/2iY = 4.0296, L/a = 220.2705, and = tt/2. The same values 
of the parameters uj p /uj h , uj p cl/c, and e as in Fig. [5] 

In the case of oblique incidence where |0| ^ 7r/2, the scat- 
tering demonstrates other interesting features, in addition to 
those discussed above. Figures [6|a) and [6jb) show the field 
distributions for the incidence angles = 3tt/7 and 6 = tt/7, 
respectively, if uj = uj i ^ uj\. The arrows in the figures 
shows the energy flow directions in the incident wave and 
in the reflected and transmitted far-zone fields. It is found 
that for 6 = 3tt/7 [see Fig.|6ja)], the reflected and transmit- 
ted rays demonstrate the behavior which is essentially differ- 
ent from the more habitual situation depicted in Fig. [6jb) for 
= 7r/7. Note that Fig.[6ja) resembles the pattern observed 
during the resonance excitation of a negative-order spatial har- 
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